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Abstract 



I report on the observation of the production of strings (dischnation hues and loops) 
via the Kibble mechanism of domain (bubble) formation in the isotropic to nematic phase 
^ • transition of a sample of uniaxial nematic liquid crystal. The probablity of string formation 

per bubble is measured to be 0.33 ± 0.01. This is in good agreement with the theoretical 
^ . value I/tt expected in two dimensions for the order parameter space S'^/Z2 of a simple 

d • uniaxial nematic liquid crystal. 
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The appearance of topological defects in phase transitions associated with sponta- 
neous symmetry breaking is common in quantum field theories, particle physics models, 
condensed matter physics and cosmology. Given a phase transition at a scale /x involving 
the reduction of explicit symmetry from a symmetry group G to a symmetry group H 

G^H, (1) 

the manifold of equivalent vacua or ground-states is the coset space A4 = G/H. Topological 
defects arise if A4 has non-trivial topology. In 3 dimensions the possibilities are classified 
the homotopy groups ito{M.) through ^^{M.). Non- vanishing ito{A4) occurs if Ai is not 
connected and leads to two-dimensional defects called domain walls. Non- vanishing 7ri(A^) 
occurs if Ai is not simply connected and leads to one-dimensional string or line defects. 
Non- vanishing n2{A4) leads to zero-dimensional point defects such as monopoles. Finally 
non- vanishing 7r3(7Vl) gives rise to defects which have been generically called texture. In 
the quantum field theory case an example is the Skyrmeon. 

Nematic Liquid Crystals (NLCs) form a class of matter with a rich variety of topo- 
logical defects. While the classification and static properties of these defects have been 
thoroughly studied [|l|-^ the mechanism of their formation and their dynamics is relatively 
unexplored in the laboratory. Recently the evolution in three dimensions of the dense 
tangle of string defects (disclinations) produced by a pressure-induced isotropic to nematic 
phase transition of a sample of NLC ("K15") was studied experimentally [P,0]- It was 
shown that the string length density (total length in string per volume) decays linearly in 
time and that the coarsening of the string density obeys the scaling hypothesis p{t) ~ t4w 
where the length scale ^(t) is assumed to be both the mean radius of curvature of strings 
and the mean separation between strings and grows in time as f^ with u = 0.51 ± 0.02. 
It was also shown that a string loop of radius R{t) shrinks in time as R{t) oc {to — t)", 
with the exponent a being 0.5 ± 0.03. This agrees with the exponent ^ obtained from an 
analysis of the nematodynamic equations of a disclination line moving through a nematic 
medium with a constant velocity v [|l|,|^. 

An important motivation for this recent work was to verify several of the string in- 
teractions postulated in the study of cosmic strings and fundamental strings. In this talk 
I report on some new work described in 0. We have observed the formation of string 
defects in the supercooling of the NLC i^l5. Preceeding the production of strings we see 
the formation of uncorrelated domains (bubbles) as described theoretically by Kibble 
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some years ago in his pioneering analysis of the production of defects (domain waUs, cosmic 
strings and monopoles) in the coohng of the early Universe. In our experiment disclina- 
tions are analogous to global cosmic strings rather than local cosmic strings as there are 
no analogues of gauge fields in simple uniaxial NLCs. The idea of testing cosmological 
theories of string formation in condensed matter systems was in fact proposed by Zurek 
Tod , who suggested a cryogenic experiment to detect superfluid helium flow arising from 



the random formation of uncorrelated domains produced by a rapid pressure quench in an 
annular geometry. 

NLCs are thermotropic materials typically consisting of rod-like molecules roughly 20 A 
in length and 5A in width. They often have two benzene rings associated with different 
functions and short flexible exterior chains. They exhibit a high-temperature isotropic 
phase characterized by disorder with respect to translational and rotational symmetries 
and a low-temperature nematic phase in which there is long-range orientational order but 
translational disorder. The nematic phase is thus liquid with respect to translations and 
crystalline with respect to orientations. Configurations of NLCs are specified by a director 
field n{x). This is a three-dimensional directionless unit vector giving at each point x the 
local orientation of the rod molecules. An order parameter for the first-order isotropic to 
nematic phase transition is ( ^^°^^ ) where 6 is the angle of the director field to some 
arbitrary fixed axis. This order parameter vanishes in the isotropic phase and is non-zero in 
the nematic phase. For K15 the isotropic to nematic transition occurs at Tnj = 35.3°C at 
atmospheric pressure. The easily accessible temperatures (typically 0° — 200°C) of nematic 
phase transitions is one reason these systems are simple to work with in the laboratory 
(especially compared to superfiuid helium). 

The effective energy density describing NLCs is similar to that of an 0(3) non-linear 
sigma model except that the NLC system is invariant only under simultaneous spatial and 
internal rotations rather than independent spatial and internal rotations [^. This means 
that the effective energy density for NLCs consists of three terms of different symmetry 
associated with splay, twist and bend deformations respectively P], as well as a surface 
term. Corresponding to each term is a potentially different elastic (Frank) constant. The 
non-linear sigma model field theory is recovered only in the isotropic limit in which all 
three elastic constants are equal (the one-constant approximation). 

NLCs exhibit three classes of topological defects. The space of ground-state nematic 
configurations A4 (ground-state manifold) is given by SO{3)/Doo, where SO{3) is the 
symmetry group of the disordered (isotropic) phase and Doo is the unbroken subgroup 



consisting of rotations about the molecular axis and 180° rotations about axes perpen- 
dicular to the molecular axis. This space is isomorphic to S'^/Z2, the two-sphere with 
antipodal points identified, which in turn is isomorphic to the projective plane RP . The 
existence of topological defects is determined by the homotopy groups of the ground-state 
manifold. The non- vanishing homotopy groups of Ai are ni{Ai) = Z2, n2{Ai) = Z and 
T^si-M.) = Z. Thus there are Z2 disclinations (string or line defects), integrally-charged 
monopoles (point defects) and integrally charged texture defects (whose field theory ana- 
logues are Skyrmeons) associated with non-trivial elements of the first, second and third 
homotopy groups respectively 0. We will be concerned here only with the tti defects. 
Note that non-trivial strings arise from the lack of polarity of liquid-crystal molecules 
which gives rise to the director field identification n = — n. This is encoded in the Z2 
isotropy group of A^ = 5'^/Z2. A path in M. which winds from n to — n as one goes 
around a closed loop in space will not be continuously deformable to a point and hence 
will represent a non-trivial disclination (string) defect. 

To see more clearly the director field configurations associated with disclinations we 
now give a very brief description of them. The basic string defect along the z-axis is 
characterized by the planar configuration of the director field 

Ux — cos(/), Uy — sin(^ and Uz — 0, (2) 

where we have taken the director to lie in the x — y plane. For an isotropic NLC (all elastic 
constants equal) the equations of motion minimizing the free energy are 

V^(i) = 0. (3) 



Solutions independent of the radial coordinate r = yx^ + y^ are given by 

(j) = sa + c, (4) 

where a = tan~^(y/x) and c is a constant between and it. In the core of the string 
(of typical radius lOOA) the director changes discontinuously and the most reasonable 
assumption is that the system is in the isotropic phase |^. The parameter s labels the 
strength of the disclination. As the core of the string is encircled the director rotates by an 
angle 27rs. The basic disclinations are the type-±i disclinations. They can be continuously 
rotated into each other and so are topologically equivalent. Integral strength disclinations 
are continuously deformable to trivial configurations without a singular core by escape into 



the third dimension and hence are topologically trivial although their dynamics is still of 
considerable interest. The energy density per unit length of a disclination is proportional to 
the square of its strength and hence the strength ±i disclinations dominate the statistical 
mechanics of these systems. 

The formation of topological defects in a phase transition as described by Kibble 
in the context of the early Universe, arises due to the existence of uncorrelated regions of 
space (domains). As these domains come into contact with each other, the variation of the 
order parameter field from one domain to another may be such that topological defects 
get trapped at the junctions of these domains. The probability of defect formation can be 
calculated by assuming that the order parameter field varies randomly from one domain to 
another. Let us illustrate this by considering a simpler case of 2 spatial dimensions when 
the ground state manifold is a circle S^. The value of the order parameter field in the 
ground state is then given by an angle 6 between and 27r. The elementary strings will 
arise when 6 winds by 27r around a closed path in space. Consider now the situation when 
three domains meet at a point and further that 6 takes random values (say ^i, 62 and ^3) 
in these three domains. The value of 6 in between any two domains will be governed by 
energetic considerations and will be such that it amounts to smallest variation of 6 from 
one domain to another. With this, the string formation at the junction will correspond to 
the case when one of the angles (say 6*3) is such that 6*1 + tt < ^3 < ^2 + tt. [Assuming 
that 62 + 1^ > di + 'n:, otherwise the order of 61 and 62 should be reversed. All angles here 
are taken mod(27r).] The maximum and minimum values of this angular span available 
for ^3 for string formation are tt and with an average of ix /2. As 6^ can assume any 
value between and 27r, the probability that a string will form at the junction is 1/4. 
The calculation of probability of formation of other defects (and for other ground state 
manifolds) can be carried out along similar lines. For NLCs the ground state manifold is 
RP"^ and the probability of string formation in 2 spatial dimensions can be calculated to 
be I/tt (see ^). 

It is known that for NLCs, the isotropic-nematic phase transition is of first order [|T2| 
and proceeds through the proccess of bubble nucleation. As a sample of NLC supercools 
below T/v/ it becomes favorable for domains (bubbles) of true ground state (nematic phase) 
to form in the false ground state medium (isotropic phase). Very small bubbles have 
too much surface energy and collapse again. Bubbles beyond a certain size (the critical 
bubbles) gain enough bulk energy to overcome the surface energy and begin to grow (see 
||13|| for a discussion of vortex formation by bubbles). The director field will randomly 



take values in Ai in different bubbles and will be roughly uniform inside a given bubble 
(for energetic reasons). As the bubbles collide, the actual correlation length will depend 
on the bubble size as well as on the rate of bubble collision and when bubble collisions 
are frequent, each bubble will represent an uncorrelated domain as used in the Kibble 
mechanism. The correlation length here is similar to the horizon in the cosmological setting 
in which different domains are causally disconnected and hence cannot communicate with 
each other. As the system cools the initially chaotic variation of the director field will 
smoothen as various neighboring domains coalesce and attempt to lower gradient energies 
by aligning their associated director fields. As described above, however, there may be 
topological obstructions to neighboring domains aligning completely. If there are non- 
trivial paths in Ai traced out by neighboring domains then tti ( string or disclination) 
defects are trapped in the system somewhat like fiux tubes in Type-II superconductors. 
This is the process we have observed in the experiment described below. One of the key 
quantities of interest in this mechanism of string production via bubbles is the probability 
of string formation per bubble. We have determined the typical domain (bubble) size 
and net string length produced thereby experimentally determining the value of the string 
formation probability. 

For our observations we used the NLC K15 (4-cyano-4'-n-pentylbiphenyl) manufac- 
tured by BDH Chemicals, Ltd.. We employed an Olympus, Inc. model BH phase contrast 
microscope with a lOX objective. The microscope was equipped with a monochrome tele- 
vision camera and a standard videocassette recorder. 

We tried two different approaches to observing the isotropic - nematic phase transition. 
In the first a freely suspended film of K15 was created by dipping a wire loop 3 mm in 
diameter in a sample of -ff 15. This wire loop was soldered to an 180 resistor which acted 
as a heater. The loop and heater assembly were fastened to an aluminum slide with a 
keyhole machined to accomodate them. With this apparatus we first heated the wire 
loop until the film entered the isotropic phase; in this phase the microscope image was 
completely structureless. We then switched off the electrical heater and allowed the film 
to cool. A similar approach was used previously by Chuang et al for the observation of 
the evolution of string defects, which formed very rapidly as the film entered its nematic 
phase. Our results in this regime largely reproduced this previous work. 

This apparatus proved to be unsuitable for observing homogeneous bubble formation 
and the subsequent generation of strings. A circular front of nematic phase formed at the 
boundary of the film and propagated through the central region of the film before we could 



finish observing the evolution of the homogeneously nucleated bubbles. We tried to reduce 
this effect by slowing down the cooling rate by slowly reducing the current through the 
resistor and by shining an auxiliary illuminator during the measurements. This improved 
the production of bubbles but did not eliminate the problem of the circular front. We 
therefore tried observing the phase transitions by using a different method by using a 
thin drop of K15 placed on an untreated, clean microscope slide. For this arrangement 
heating was done using an illuminator; we were able to obtain satisfactory images of bubble 
formation and evolution. 

One set of such images is reproduced in Fig. 1. Fig. 1(a) shows the numerous small 
isolated bubbles of nematic phase which form first. Fig. 1(b) and 1(c) show the images 
at short intervals later, as the nematic bubbles increase in size. This occurs both by 
natural growth and also by coalescence. In the next stage the organization of the NLC 
into bubbles is replaced by an image of a homogeneous medium with entangled strings 
(Fig. 1(d)), which further evolve by straightening, shrinking, and excising small loops of 
closed string (Fig. 1(e)). The string dynamics has been very well described previously 
1^,0. The film becomes increasingly homogeneous in appearance on longer time scales as 
the string pattern coarsens. 

We analyzed series of images such as Fig. 1 to estimate the string to bubble ratio. In 
counting the number of bubbles we have to decide which bubbles are relevant for string 
formation. For example between Fig. 1(b) and 1(c) small bubbles seem to just get absorbed 
by bigger bubbles and we do not see them contributing to any string formation. This 



observation is consistent with the results of the numerical simulation in |]T^ where the 
collision of two large (critical) bubbles with one small (subcritical) bubble was considered 
and it was found that even if a string forms in such a collision, it gets pushed out of the 
smaller bubble due to large momentum of the fast expanding walls of the large bubbles. 
Even in the collision where all the bubbles involved are very small, we have not observed any 
string formation. This could again be because a string will get trapped at the collision point 
only when the bubble walls have sufficient momentum which may not happen if bubbles 
are very small. Even if some strings do form in the collisions of very tiny bubbles, these 
strings will be very small (for isolated collisions) and will disappear quickly by shrinking. 
As we count the strings at a little later stage anyway (when we can clearly distinguish 
them), these small strings would have disappeared by that time and hence no error will be 
caused by the neglect of very small bubbles. 
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Guided by these considerations, we made a judgement of a typical size for the bubble 
which seems to be large enough to contribute to string formation and then consider coa- 
lescence and growth of such bubbles. Bubbles r times the average bubble size are assumed 
to have arisen from the coalescence of r average bubbles (as can be verified by looking at 
the earlier stages of the bubble formation) and are thus counted r times. 

The next step is to look at a photograph corresponding to the time at which these A^ 
bubbles have expanded to fill up the entire region of the photograph (of area A). As the 
bubbles coalesce they do not remain spherical but actually squeeze into gaps. The final 
image looks roughly like a rectangular region filled by A^ squares (with typical side of a 



square d roughly equal to ^/AJN). We estimated a count of about 55 "effective" bubbles 
in Fig. 1(c); the typical length scale was about 80 /um. 

Our estimate of the probability that each of these bubbles will lead to a string upon 
final coalescence was made as follows. Imagine that the rectangular image above is a lattice 
of area A with m rows and n columns such that A^ = mn. As strings form because of 
mismatches of the director at bubble interfaces the maximum possible string length L^ax 
is given by the total sum of the lengths of all the links on the lattice which is 2{NAY''^. 
For Fig. 1(c) we estimate Lmax = 8.7 mm. We then examine a photograph of the earliest 
possible stage when one can discern strings clearly and measure the total length of string 
Ls'i in Fig. 1(d) we obtained 2.7 mm. The probability of string formation is then taken to 
be P = Ls/Lmax = 0.31. 

We repeated this analysis on several sequences, obtaining P = 0.33±0.01 for the string 
formed per bubble. Systematic errors are not negligible; we believe that this measurement 
provides a lower bound on the true value of P for two reasons. First the total string length 
will have decreased from its initial value by the time we are able to clearly resolve and 
measure Lg. Secondly strings which wander significantly in the direction normal to the 
plane of focus (but are still formed due to top bubble layers which we observe) will have 
their length underestimated. The actual underestimation may not be too bad since string 
wiggling will be dampened in the time between their formation and our measurement of 
their length. Further there is a selection effect in the experiment since we are able to 
see strings clearly only when they do not deviate significantly from the plane of focus. 
Thus we are effectively counting planar strings. This will suppress the effects of the three- 
dimensionality of the bubble layers. In other words, the bubbles we see are only the top few 
layers; similarly we see the strings at the time they are roughly confined to these top layers. 
The numbers we obtain, therefore, correspond to an effectively two-dimensional situation. 



Note that for strictly 2 dimensional situation, the strings will form perpendicular to the 
plane formed by the bubble layer. In the present case, however, strings will bend around in 
the plane as they come towards the top (where the director will eventually become normal 
to the surface this being the preferred boundary condition for NLC-air interface) as well 
as when they squeeze through the bubble layers underneath them. 

We note here that the string to bubble ratio we have measured is in good agreement 
with the theoretical value I/tt ^ 0.32 obtained by Vachaspati [^ for the ground-state 
manifold S"^ /7j2 in two dimensions, especially considering the above arguments that our 
estimates for P correspond to an effectively 2 dimensional situation. 

We now turn to another interesting phenomenon we are able to observe in this exper- 
iment. We first cool our sample through the nematic phase transition so that disclinations 
are produced and then superheat from the nematic phase to the isotropic phase. Spherical 
isotropic bubbles are produced in the nematic medium. These have a somewhat different 
more transparent appearance than the nematic bubbles in an isotropic medium because of 
the different optical properties of the two phases. In this way it is quite easy to distinguish 
the different types of bubble. We observe very clearly in this process the heterogeneous 
nucleation of isotropic bubbles along disclination lines like beads on a necklace (see Fig. 2). 
This nucleation process is of great technological and theoretical interest |jl5,|15[ . We feel 
that NLCs are a good testing ground to check experimentally many of the theoretical pre- 
dictions concerning the mechanism of the nucleation process such as the size distribution 
of bubbles. 

Finally we point out that the implications of the work described here extend beyond 
the dynamics of the formation of disclinations in nematic liquid cystals because of the close 
analogies with the formation of other one-dimensional topological defects such as screw 
dislocations in crystals and vortices in superfluid helium P,p!0[|. 

First of all I wish to thank my collaborators L. Chandar, Eric Schiff and Ajit Srivas- 
tava. Thanks are due to Bernard Yurke for his valuable advice during his visit to Syracuse 
University. I would also like to express very warm thanks to Edward Lipson and his bio- 
physics group at Syracuse University for their generosity, both for sharing their microscope 
apparatus and for instruction in its use. Without their support the research reported here 
would not have been possible. This research was supported by the Outstanding Junior 
Investigator Grant DOE DE-FG02-85ER40231. Finally I am indebted to Norma Sanchez 
for her kind invitation to lecture at Erice. 
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Figure Captions 

Fig. 1. A series of five images showing the isotropic to nematic phase transition in a drop 
of K15 on an untreated microscope slide; the images show a width of 0.7 mm of 
the drop. Note the stages of bubble nucleation and growth, bubble coalescence 
and string formation, and string coarsening. The delay times for each image 
(referred to the first frame showing discernible bubbles) are (a) 2 s, (b) 3 s, (c) 5 
s, (d) 11 s and (e) 23 s. The scale of the pictures is given in (f). 

Fig. 2. This series of three pictures shows the heterogeneous nucleation of isotropic bub- 
bles along disclination lines. The delay times for each image, referred to the 
picture in (a), are (b) 2 s and (c) 4 s. The scale is given in (d). 
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